Neutron drip line and the equation of state of nuclear matter 
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We investigate how the neutron drip line is related to the density dependence of the symmetry 
energy, by using a macroscopic nuclear model that allows us to calculate nuclear masses in a way 
dependent on the equation of state of asymmetric nuclear matter. The neutron drip line obtained 
from these masses is shown to appreciably shift to a neutron-rich side in a nuclear chart as the 
density derivative of the symmetry energy increases. Such shift is clearly seen for light nuclei, a 
feature coming mainly from the surface property of neutron-rich nuclei. 
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Thanks to recent developments of radioactive ion beam 
facilities, one might be able to experimentally probe the 
stability of atomic nuclei against neutron or proton drip 
in a nuclear chart ranging from light to superheavy nu- 
clides. The key quantities to study the neutron (proton) 
drip line are the one- and two-neutron (proton) separa- 
tion energies, S^tp) and S2n{2p): which correspond to an 
energy required to remove one and two neutrons (pro- 
tons) from a nucleus in the ground state, respectively. 
Experimentally, the neutron drip line is marginally acces- 
sible only for light nuclei [l|. Even beyond the neutron 
drip line, however, nuclei can be present in dense neutral 
matter. Nuclei in the crust of neutron stars are a typical 
example and are considered to control the thermal and 
electric transport properties of matter in the crust as well 
as the dynamics of superfluid neutron vortices, which are 
relevant to the observed thermal and rotational evolution 
of neutron stars 2| . We remark that the size and shape 
of nuclei are shown to be controlled by the equation of 
state (EOS) of asymmetric nuclear matter through the 
density dependence of the symmetry energy [3|] . In this 
Letter, we will investigate how the density dependence 
of the symmetry energy in turn affects the prediction of 
the neutron drip line. 

Theoretically, a Weizsacker-Bethe mass formula, which 
is based on a view of nuclei as incompressible spherical 
liquid drops of uniform density uq, provides a standard 
behavior of the neutron drip line. In this formula, the 
nuclear binding energy Eb is written as function of mass 
number A and charge number Z (or neutron number N) 
in the form 
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where E^oi = a^oiA is the volume energy, Egym 
asvm[(^ — Z)/A]'^A is the symmetry energy, i^surf 



Osurf^^^"^ is the surface energy, and i^coui = 
acouiZ^/A^/^ is the Coulomb energy. Then, the one- 



neutron separation energy can be evaluated as 
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where x — Z/A is the proton fraction. The condition 
Sn = gives a smoothed behavior of the neutron drip 
line. The derived drip line is close to x = 0.3, which 
is basically controlled by the competition between the 
volume and symmetry energy terms, except in the light 
region of the nuclear chart. This may be a good start- 
ing point, but one needs to go beyond Eq. ^ by taking 
into account nonnegligible deviation of the nuclear den- 
sity from hq. 

Recently, the density dependence of the symmetry en- 
ergy attracts much attention because it is relevant to the 
isospin dependence of nuclear masses (e.g., Refs. [4, 5)]) 
and radii (e.g., Refs. [6|-l8|), dipole resonances (e.g., Refs. 
[Sl llOj). and heavy- ion collisions involving neutron-rich 
nuclei (e.g., Refs. |lll4l3| ). In predicting the neutron 
drip line, uncertainties in the density dependence of the 
symmetry energy need to be taken seriously. The impor- 
tant parameter characterizing the density dependence of 
the symmetry energy is a density symmetry coefficient L, 
which is defined as L = 3no{dS / dn)n=no with the sym- 
metry energy S{n) dependent on the density n of bulk 
nuclear matter. Masses of extremely neutron-rich nuclei 
were calculated from a macroscopic nuclear model and 
shown to have an appreciable dependence on L, which 
can be understood from the density and isospin depen- 
dence of the surface tension [5| . Here we address how this 
dependence affects the neutron drip line on the nuclear 
chart. 

We begin with a macroscopic model of nuclei [6| , which 
was constructed in such a way as to reproduce the known 
global properties of stable nuclei and can be used for 
describing the masses and radii of unstable nuclei in a 



manner that is dependent on the EOS of nuclear matter. 
This model can be summarized as follows: 
(i) We set the bulk energy per nucleon as 
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are the potential energy densities for symmetric nuclear 
matter and pure neutron matter, n„ and Up are the 
neutron and proton number densities, n == n„ + Tip, 
ct = {rin — np)/n is the neutron excess, and m„ is the 
neutron mass. A set of expressions ©-([S]) is one of the 
simplest that reduces to a usual expansion 14 1 
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in the limit of n — >■ no and a ^ Q. Here wq and Kq 
are the saturation energy and the incompressibility of 
symmetric nuclear matter, and Sq = S{n = uq). In 
the incompressible limit, wo and S'o correspond to a vol 
and asym in the mass formula ([T]), respectively. We fix 
63, which controls the EOS of matter for large neutron 
excess and high density, at 1.58632 fm'^. This value was 
obtained by one of the authors 15[ in such a way as to 
reproduce the neutron matter energy of Friedman and 
Pandharipande [l6[. Change in this parameter would 
make no significant difference in the determination of the 
other parameters and the final results for nuclear masses, 
(ii) We write down the total energy of a nucleus of mass 
number A and charge number Z as a function of the 
density distributions n„(r) and n-p(r) in the form 



E = Eb + Eg+Ec + Nninc^ + ZnipC^, 



where 
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is the bulk energy. 
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is the gradient energy with adjustable constant Fq, 
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FIG. 1: (Color online) The neutron and proton drip lines 
obtained from the EOS models C and G by using the macro- 
scopic nuclear model and from a contemporary mass formula 
[17| . The regions filled with squares correspond to empirically 
known nuclides [13, fTyl. 



is the Coulomb energy, and rup is the proton mass. 

(iii) For simplicity we use the following parametrization 

for the nucleon distributions ni{r) {i — n,p): 
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where r is the distance from the center of the nucleus. 
This parametrization allows for the central density, half- 
density radius, and surface diffuseness for neutrons and 
protons separately. 

(iv) In order to construct the nuclear model in such a 
way as to reproduce empirical masses and radii of sta- 
ble nuclei, we first extremize the binding energy with 
respect to the particle distributions for fixed A, five EOS 
parameters, and Fq. Next, for various sets of the incom- 
pressibility and the density symmetry coefficient, we ob- 
tained the remaining three EOS parameters and the gra- 
dient coefficient by fitting the calculated optimal values 
of charge number, mass excess, root-mean-square (rms) 
charge radius to empirical data for stable nuclei on the 
smoothed /3 stability line |13| . In the range of the param- 
eters < i < 160 MeV and 180 MeV < Kq < 360 MeV, 
as long as KoSo/3noL > 200 MeV fm^, we obtained a 
reasonable fitting to such data. As a result of this fit- 
ting, the parameters uq, wq, Sq, and Fq are constrained 
as no = 0.14-0.17 fm-^, wq == -16 ± 1 MeV, S'o = 25-40 
MeV, and i^o = 66 ± 6 McV fm^. The fitting gives rise 
to a relation nearly independent of Kq, 



Sq^B + cl, 



(12) 



where B « 28 MeV and C « 0.075. 

We proceed to obtain the neutron and proton drip 
lines from the macroscopic nuclear model. For various 
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. 2: (Color online) N — Z obtained for nuclides at neutron 
proton drip from the EOS models C and G. 
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FIG. 3: (Color online) Differences in the neutron number A'^ 
obtained for nuclides at neutron and proton drip between the 
calculations from the EOS models C and G. 



sets of L and Kq, we first evaluate the binding energy 
Eb of nuclei in the ground state by minimizing the en- 
ergy d?]) for fixed N and Z. We then draw the neu- 
tron (proton) drip line by identifying nuclides at neu- 
tron (proton) drip with those neighboring to nuclides for 
which Sn = EBiZ,N)-EB{Z,N-l) {Sp - EBiZ,N)~ 
Eb{Z - l,N)) and S2n = Eb{Z,N) - Eb{Z,N - 2) 
{S2p = Eb{Z, N) - Eb{Z - 2, N)) are positive and be- 
yond which at least one of them is negative. The results 
obtained from the two extreme EOS models denoted as 
EOS C (L = 146 MeV and K^ = 360 MeV) and EOS 
G (L = 5.7 MeV and K^ = 180 MeV) are shown in Fig. 
1 , together with the empirically known nuclides 18J, |l9[ 
and the prediction from a contemporary mass formula 
17| . We remark that inclusion of the condition for S2n 
and S2p in addition to Sn and Sp in drawing the drip lines 
makes only a little difference in the case of the present 
model calculations, while being essential in the case of the 
prediction from the mass formula because of the Wigner, 
shell, and even-odd terms included therein. We remark 
that the rms deviations of the calculated masses from the 
measured values [10| are about 3 MeV, which is of order 
the deviations obtained from a Weizsacker-Bethe mass 
formula. 

We find from Fig. 1 that the obtained neutron drip 
lines show an appreciable L dependence, while the pro- 
ton ones do not. This is reasonable because nuclides at 
neutron (proton) drip are far away from (near) N = Z 
(see Fig. 2). The neutron drip line does shift to a neutron- 
rich side as L increases, a feature that will be discussed 
later in terms of a compressible liquid-drop model. We 
remark that the EOS dependence of the obtained drip 
lines comes predominantly from L because of negligible 
Kq dependence of the calculated masses [5| . 

In order to see the L dependence more clearly, we plot 
in Fig. 3 the difference in the neutron number of nuclides 
at neutron and proton drip between the calculations from 
the EOS models C and G. The difference shows only a 
weak dependence on Z both in the case of neutron and 
proton drip. This indicates that the L dependence can 
be seen more clearly for lighter nuclei. In fact, the cor- 
responding proton fraction of nuclides at neutron drip 
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FIG. 4: (Color online) Proton fraction for nuclides at neutron 
drip obtained from the EOS models C and G. 



shows a stronger dependence on L for lighter nuclei, as 
shown in Fig. 4. This is advantageous because heavier ra- 
dioisotopes are more difficult to produce in experiments. 
The L dependence of the neutron drip line as obtained 
above can be understood within the framework of a com- 
pressible liquid-drop model in which nuclei in equilibrium 
are allowed to have a density different from the saturation 
density no of symmetric nuclear matter. By following a 
line of argument of Ref. [5| , we first add the surface sym- 
metry term, assymA'^^^liN — Z) /A]'^ , to the mass formula 
p]) based on an incompressible liquid-drop model. This 
surface symmetry term gives rise to additional contribu- 
tion, 
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to the neutron separation energy ([2]). Next, we consider 
the density-dependent surface tension [20|, 
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are the density and neutron ex- 
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the surface symmetry energy coefficient, and x — 
(no/o'o)9cr/9nin|ni„=no,Qi„=o- By taking a limit of van- 
ishing compressibility, one obtains AiraoR^ = Osurf^^^"^ 
and AiraQCsymR^ = — Ossym^^ with the liquid-drop ra- 
dius R. Typically, fitting to the empirical mass data 



cess inside a liquid drop, aa = cr(no,0), v_/sym 



yields ctq « 1 MeV fm~^ and Cgym = 1.5-2.5. As we 
shall see below, nonvanishing compressibility effectively 
introduces the L dependence into the surface symmetry 
term through the parameter x characterizing the den- 
sity dependence of the surface tension. The value of x is 
poorly known, but likely to be positive [5j . For example, 
X — 4/3 for the Fermi gas model. 

If one ignores Coulomb and surface corrections, the 
equilibrium density and energy per nucleon of a liquid- 
drop. Us and Ws, can be evaluated from Eq. ([6]) as 



and 
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Strictly speaking, expressions ([T5|) and ([TB| are applica- 
ble only for nearly symmetric nuclear matter. We nev- 
ertheless use these expressions for the purpose of char- 
acterizing the liquid-drop properties because the typical 
value of a along the neutron drip line is of order 0.35-0.4, 
considerably smaller than unity. By substituting Us , Eq. 
([T6| . into nin in Eq. ([T4| . we thus obtain 
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This result allows one to identify assym^^^^ with 
— 47r(To-R^(C'sym + 3Lx/Ko) and hence to conclude that 
with increasing L, Sn increases through the surface sym- 
metry contribution (|13l) for x < 1/2. 

This conclusion is consistent with the L dependence of 
the neutron drip line shown in Fig. 1 because any positive 
corrections to Sn act to enhance the stability of nuclei 
against neutron emission. It is important to note that 
the bulk symmetry term gives rise to a negative contri- 
bution to Sn for a; < 1/2, as shown in Eq. ([2]), and that 
the parameter asym corresponds to the symmetry energy 
coefficient ^o, which in turn is related to L by the relation 
P^ obtained from fitting to empirical masses and charge 
radii of stable nuclei. Since one obtains a larger Osym for 
larger L, the effect of flgym tends to decrease S'„ with L 
and hence to facilitate neutron drip. However, this effect 
is relatively small compared with the above-mentioned 
effect of flssym. This is consistent with the fact that the 
L dependence is clearer for lighter nuclei. 

In summary, we have investigated the influence of the 
density dependence of the symmetry energy on the drip 
lines by using a macroscopic nuclear model that depends 
explicitly on the EOS of nuclear matter. We find that an 
L dependence appears appreciably in the neutron drip 
line and it is clearer for lighter nuclei, a feature com- 
ing mainly from the surface property through the den- 
sity and neutron excess dependence of the surface ten- 
sion. We note that our calculations do not include even- 
odd or shell corrections. The even-odd corrections would 



play an important role in Sn- In fact, the magnitude of 
even-odd staggering in Sn could be comparable to that of 
change in Sn due to uncertainties in L. This implies that 
some kind of smoothing would be required in deriving 
information about the EOS from future empirical data 
on neutron drip. Shell corrections would further compli- 
cate such derivation, but are intriguing in the context of 
magicity of nuclei close to the neutron drip line 2l| . 
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